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ABSTRACT 

We present a new formulation of Lagrangian perturbation theory which allows accurate pre- 
dictions of the real- and redshift-space correlation functions of the mass field and dark matter 
halos. Our formulation involves a non-perturbative resummation of Lagrangian p erturbation 
theory and indeed can be viewed as a partial resummation of the formalism of Matsubara 
(2008a b) in which we keep exponentiated all of the terms which tend to a constant at large 
separation. One of the key features of our method is that we naturally recover the Zel'dovich 
approximation as the lowest order of our expansion for the matter correlation function. We 
compare our results against a suite of N-body simulations and obtain good agreement for the 
correlation functions in real-space and for the monopole correlation function in redshift space. 
The agreement becomes worse for higher multipole moments of the redshift-space, halo cor- 
relation function. Our formalism naturally includes non-linear bias and explains the strong 
bias-dependence of the multipole moments of the redshift-space correlation function seen in 
N-body simulations. 

Key words: gravitation; galaxies: haloes; galaxies: statistics; cosmological parameters; large- 
scale structure of Universe 



1 INTRODUCTION 

The observed large-scale structure (LSS) of the universe is a pil- 
lar of modern observational cosmology, providing a window into 
the primordial fluctuations, expansion history, and growth rate of 
perturbations, as well as allowing tests of the theory of gravity 
on the largest accessible scales. The two-point correlation func- 
tion (or its Fourier transform) is a useful and relatively simple 
compression of the cosmological information of interest. How- 
ever, the interpretation of LSS statistics is hampered by two pri- 
mary uncertainties: LSS tracers (e.g., galaxies, Ly-a forest, 21 
cm) are biased relative to the underlying matter density field 
and are observed in redshift space. On very large scales these 
two effects are simpl e linear transformations of the und erlying 
matter density field jKaiserl 1 19871 ; lEfstathiou et al.l Il988h . while 
on Mpc scales, the dynamics are highly non-linear, and N-body 
simulations seem to be required for quantitative accuracy. How- 
ever, on intermediate or quasi-linear scales there is hope that 
observable quantities for biased tracers may be accurately mod- 
eled semi-analyt ically by extending perturbation theory beyond 
linear order (see.lTaruva. Nishimichi & Saitdl201Ch iReid I & White! 



201 1 ; Okamura, Tanrya, & Matsubara 201 1 ; Elia et al. 201 1 ; 



H2. 



Crocce, Scoccimarro, & Bernardeau 2012; Vlah et al. 2012, for re- 
cent work in this direction). iMatsubar j 12008 a b) introduced a new 
perturbative scheme (which we shall refer to as Lagrangian Re- 
summation Theory; LRT) which addresses both non-linear bias- 
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ing and redshift space distortions in a unified framework based 
on Lagrangian perturbation theory and that substantially improves 
upon standard perturbation theory for the description of both 
matter and dark matter halo clustering in the quasilinear regime 
( Padmanabhan & Whitell2009l: iNoh. White & Padmanabhad l2009t 
IReid & Whitel201 it I Sato & Matsubarall201 ll : lRampj2012h . In this 
pa per propose a new r esummation scheme which extends the work 
of Matsubara] J2008dlb1) and results in a more accurate expression 
for the two-point correlation function in both real- and redshift- 
space for both matter and dark matter halos. 



There are several advantages to adopting a Lagrangian de- 
scription of the L SS. The well-known Zel'dovich approximation 
dZel'dovichlll970l) provides a one time-step, reasonably accurate 
approximation to the non-linear density field by displacing La- 
grangian particles by t he linear theory displacemen t field (for a 
recent examination, see lTassev & Zaldarriag a 2012a. b). A distinc- 
tive feature of the present work is that we recover the Zel'dovich 
result exactly as a limit of our expression for the correlation 
function of matter (which is not the case in a similar study by 
iMcCullagh & Szalavl 120121) . The clustering of dark matter halos 
is of greater interest for the interpretation of galaxy redshift sur- 
veys, since modem galaxy formation models assume that galaxies 
form and reside in the gravitational potential wells of dark mat- 
ter halos. Again, the LRT approach is advantageous since a local 
Lagrangian biasing scheme provides a better description of the bi- 
asing of dark matter halos in N-body simul a tions compared with 
local Eulerian bias dRoth & Porcianil 1201 ll ; iBaldauf et all 120121: 
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IChan. Scoccimarro. & Shethll2012h an d can be extended to in clude 
a continuous galaxy formation history dWang & Szalav 20l3). 

The outline of the paper is as follows. We begin with a brief 
review of Lagrangian perturbation theory in Section[2] In Section[3] 
we present our results within the context of the Zel'dovich approx- 
imation, in which many of the main physical effects are present 
but the algebra is simplified. In Section [4] we extend our results 
to the next higher order in perturbation theory, with most of the 
details and formulae being relegated to Appendices. In Section [5] 
we present a detailed comparison of our analytic theory with high- 
precision N-body simulations. We conclude in Section[6] 

For plots and numerical comparisons we assume a ACDM 
cosmology with Q. m = 0.274, fi A = 0.726, h = 0.7, n = 0.95, 
and o" 8 = 0.8. Our simulation data are derived from a suite 
of 20 N-body s imulations run with the TreePM code described 
in IWhitel j2002l) . Each simulation employed 1500 3 equal mass 
{m p =* 7.6 x 10'° /r'M Q ) p articles in a periodic cube of side lengt h 
1.5 /r'Gpc as described in dReid & Whitel201 lMWhite et al.l201ll) . 



2 BACKGROUND AND REVIEW 

In this section we provide a brief review of cosmologica l pertur- 
bation theory, focusing on the Lagrangian fo r mulatiorFl dBuchert 
1989 1 ; iMoutarde et aill99ll ; Iffivon et alj|l995l : iTavlor & Hamilton 
1996). This material should be sufficient to remind the reader of 
some essential terminology, and to establish our notational con ven- 
tions. Our discussion is largely drawn from lMatsubaral q2008al lbf) to 
which we refer the reader for further details. 



the Fourier transform convention 



2.1 Basic definitions 

Cosmological perturbation theory concerns itself with predicting 
the clustering properties of cosmological fluids. In the context of 
large-scale structure, (i.e., at late times when baryons and pho- 
tons have completely decoupled), the only relevant fluid is the mat- 
ter fluid, which, on all but the very smallest scales, interacts only 
through self-gravitational coupling. 

The matter fluid is idealized as a single-streaming, pressure- 
less dust, characterized at any time t by its mass density p(x, t) and 
peculiar velocity field v(x, t). Following common convention, we 
let x denote position in comoving coordinates, and t denote proper 
time for a comoving observer. The mean value p(t) = (p(x, t)) of 
the mass density decreases as the universe expands like p oc a~ 3 , 
where a{i) is the cosmic scale factor. (Here and hereafter, we drop 
the explicit time dependence in equations where all quantities are 
to be evaluated at the same time.) Deviations from homogeneity 
are expressed in terms of the density contrast S(x, t), defined by the 
relation p(x) = p[l + 6(x)]. 

The most important statistical quantities that can be formed 
from these fields are the 2-point correlation function, 

m = (6(x)6(x + r)>, (1) 

and its Fourier transform, the power spectrum P(k), defined by 

(S(k)S(k')) = (2n) 3 6 D (k + k')P(k). (2) 

Here S D denotes the 3-dimensional Dirac delta function, and we use 
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F(x) 



f- 

J (2n 



i F(k)e ik *. 



(3) 



Angle brackets around a cosmological field, e.g. (F), signify an en- 
semble average of that quantity over all possible realizations of our 
universe; in most cases of interest, ergodicity allows us to replace 
these ensemble averages with spatial averages over a sufficiently 
large cosmic volume. 



2.2 Lagrangian perturbation theory 

In the Lagrangian approach to cosmological fluid dynamics, one 
traces the trajectory of an individual fluid element through space 
and time. For a fluid element located at position q at some initial 
time to, its position at subsequent times can be written in terms of 
the Lagrangian displacement field *P, 



x(q, t) = q + ¥(9, f), 



(4) 



where *F(<jr, ?o) = 0. Every element of the fluid is uniquely labeled 
by its Lagrangian coordinate q, so that for a fixed / we may view 
the mapping q «-» x as a simple change of variable. 

The displacement field *¥(q,t) fully specifies the motion of 
the cosmological fluid. The aim of Lagrangian Perturbation Theory 
(LPT) is to find a perturbative solution for the displacement field, 



T(9, r) = *¥ m (q, + "F <2) (9, t) + >P (3) (9, t) + ■ 



(5) 



The first order solutio n is the well-known Zel'dovich approxima- 
tion dZerdovichlll970l). Explic it solutions are known up to fourth 
order dRampf & Buchertll2012l) . 
The continuity equation 



[1 + 6(x,t)]d 3 x = [l+S (q)]d 3 q 



(6) 



expresses the fact that, for a smoothly evolving fluid, an element 
d 3 q centered at q at time t is transformed into an element d 3 x cen- 
tered at x(q, t) at time t. The initial time t may be taken to be early 
enough that the initial matter fluctuations do(q) are arbitrarily small, 
so that we may formally express the Eulerian density field in terms 
of the Lagrangian displacement field as 



1 + S(x, t) : 



q6 D [x-q-'¥(q,t)]. 



(7) 



2.3 Biased tracers 

Although small, the initial density fluctuations 6o(q) provide the 
seeds for subsequent structure formation. In this paper we restrict 
ourselves to a local Lagrangian bias model, which posits that the 
locations of discrete tracers at some late time t are determined by 
the overdensities in the initial matter density field. More explicitly, 
let Sniq) denote the matter density contrast at the initial time to, 
smoothed on some scale R. (In the end the value of this smoothing 
scale will turn out to be irrelevant, but its use helps ensure that 
intermediate quantities are well-behaved.) Consider a collection of 
discrete tracers X, where X might denote galaxies of a particular 
type, or halos of a particular mass range, etc. The locations of 
these tracers at time t may be identified with particular points (9,) 
in the initial density field by inverting Eqn. [4] Our hypothesis is 
that these initial locations {q t ) are drawn from a distribution that is 
a locally biased function of the smoothed matter density field, i.e. 



Px(q) = PxF[S R (q)]. 



(8) 
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Here p x is the mean comoving number density of our tracer X 
and the function F{5) is called the Lagrangian bias function. The 
perturbations in p x (q) are 0(6 R (q)) and therefore also arbitrarily 
small. When this biasing relation is viewed in Eulerian coordinates, 
it is non-local: p x (x) is depends on the matter densi ty at points 
other than x (e.g. JCatelan et alTl 998; Matsubara 201 1). The corre- 
sponding non-local Eulerian bias terms can b e most easily seen in 
their contribution to the bispectrum of halos dBaldauf et alj|2012l ; 
Ch an. Scoccimarro. & Shethll2012h.|Matsubard J201 lh extends the 
formalism from lMatsubaraH2008al lbl) that we have adopted here 
to incl ude non-local Lagran gian bias as well, e.g., the peaks bias 
model l lBardeenetaTi ri986), but we do not explore those exten- 
sions here. 



2.4 Redshift space 

While analyzing the clustering of biased tracers is difficult enough, 
for modern surveys we must also deal with the complication of 
redshift space distortions. The position of an object, located at true 
comoving position x, will be mis-identified due to its peculiar ve- 
locity along the line-of-sight, as 



(9) 



In this work we adopt the standard "plane-parallel" or "distant- 
observer" approximation, in which the line-of-sight direction to 
each object is taken to be the fixed direction z. While this may seem 
a poor assumption for modern wide-area surveys, it has been shown 
to be sufficient within the level of current e rror bars (e.g., Figure 10 
of lSamu shia. Per cival. & Raccanellill2012h . 

In the Lagrangian approach, including redshift-space distor- 
tions requires only a simple additive offset of the displacement 
field. The peculiar velocity of a fluid element, labeled by its La- 
grangian coordinate q, is at any time given by 



v(q) = ax(q) = d9(q). 



(10) 



So in redshift space, the apparent displacement of the fluid element 
is 



z ■ *P 
1* = *P + : — z. 

H 



(11) 



To a good approximation the time dependence of the nth order term 
in Eq. © is given by T " oc D". Therefore * ( "' = nHfY"\ where 
/ = dlogD/dloga is the growth rate, often approximated as / ss 
flJJ, 6 . The mapping to redshift space may then be achieved, order- 
by-order, via the matrix 



F$ = 6,j + nfi,tj. 



(12) 



3 ZEL'DOVICH APPROXIMATION 

In this section we present a derivation of our new result in 
the simplified setting of the Zel'dovich approximation (see also 
iBond & Couchman 1988). This allows us to sketch the main idea 
of our approach while avoiding many of the complications inherent 
in perturbative calculations. Several key points regarding the form 
of the solution are made along the way. 

Our starting point is the continuity equation 



[1+S x (x,t)]d 3 x=[l+S x (q,t )] d 3 q, 



expressing the conservation of number density for the tracer X be- 
tween times to and t. Invoking the hypothesis of local Lagrangian 
biasing, the quantity on the right-hand side is 



so that 
1 + 6 x (x, t) 



l+S x (q,t Q ) = F[S R (q)], 
dx 



dq 



F[S R (q)] 



d 3 q F[6 R (q)]6 D [x-q-<¥(q,t)] 



(14) 

(15) 
(16) 



In the following, we will suppress the explicit dependence on t 
when there is no risk of ambiguity. We now replace the delta func- 
tion with its Fourier representation, and also introduce the Fourier 
transform F(A) of F(S), 

I +<Hm= ! a\,FU h 4 ( p\ | £L t .« -/-^i ,17, 



■S x (x) = J d 3 q F[S R (q)} J -^ e ' Hx ~ 



(18) 

The 2-point correlation function (; x (r) = (6 x (xi)S x (x2)) for 
the biased tracer X is then given by 

l+£,(r) = f> gl d\ 2 f f \ f \ 
J J (2tt) 3 (2tt) 3 

where 6 a = 6 R (q a ), *P„ = ¥(<?„), and r = xi - x\. By statisti- 
cal homogeneity, the expectation value above depends only on the 
difference in Lagrangian coordinates, q = q 2 - q l . The change of 
variables {q lt <jr,) — > [q,Q = (q x + q^)l2} then leads to 

■I 3 '- ... r dA, riJ-, . . 

(20) 



/, C d k ■.,„ . C dA\ dA-> - ~ 
d q i (In) 3 J ~7jt^2jt FlF ~ 



where we have defined 

K{q,k,A u A 2 ) 



i(Ai5i+A 2 d2+k 



A) >> 



(21) 



and A = *P 2 - *Pi . This expression is th e exact configuration space 
analog of Eq. (9) in lMatsubar J j2008bh . 

The cumulant expansion theorem allows us to expand the ex- 
pectation value in Eq. J2U in terms of cumulants, 



(e ) = exp 



(22) 



where (X N ) C denotes the Nth cumulant of the random variable X. 
The field S R (q) is a smoothed version of the linear density field 
S L (q), and is therefore Gaussian. Within the Zel'dovich approxima- 
tion, the displacement field 



f(q) 



J (2n) 3 k* 



Sdk), 



(23) 



is linear in S L , hence also Gaussian. Thus, in the applying the cu- 
mulant expansion theorem to Eq. J21b . only the second cumulant 
survives, 

= (A]+A 2 2 )o 2 R +A lJ k l k 1 
+ 2A l A2^R + 2(A l +A 2 )U,k„ (24) 



((A^! + A 2 S 2 + k ■ A) 2 ) c 



(13) 



where we have defined 
A,j(q) = (AiAj) c , 



Uq) = <*ife>« (25) 
Ui(q) = (6^ = {6 2 Ai) c . (26) 
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Eq. l l2 1 1 then evaluates to 



K = exp 



(27) 

The quantity a\ is simply the variance of the smoothed lin- 
ear density field, while = (<5«(<Zi)<5/;(92)) is m e corresponding 
smoothed linear correlation function. The matrix Ay may be de- 
composed as 



(2S) 



where cr 2 



■^<|*P| 2 > is the 1-D dispersion of the displacement field, 
and 7711 and n ± are the transverse and longitudinal components of the 
Lagrangian 2-point function, rjijiq) = (^(ffi)^/^))- The vector 
Uj(q) = U (q) q t is the cross-correlation between the linear density 
field and the Lagrangian displacement field. In the Zel'dovich ap- 
proximation these quantities are given by 



»7||0?) 



q±(q) 

1 

U(q) = 



2^L 



dk P L (k), 

h(kq) 



dk P L (k) 



dk P L (k) 



jo(kq) - 2 



j\(kq) 



dkkP L (k) j\(kq). 



(29) 



(30) 



(31) 



(32) 



Up to factors of 2 and /, these expressions are ide ntical to the 
Eulerian velocity c orrelators in linear theory (e.g. iFishen 1 19951 ; 
iReid&WMtekoill) . which is not surprising since v L = f¥ in the 
Zel'dovich approximation. 

3.1 Exact results for matter 

At this point we pause to consider the unbiased case, where F(S) = 
1 or F(A) = 2nS D (A). In this limit Eq. d20t reduces to 



1 +f 



\Aijkjkj 



I; 



d 3 q 



-^(r-q^A-Hr-q) 



(33) 



(34) 



(2nfl 2 \A\ m 

after carrying out the Gaussian integral over k analytically. This is 
an exact expression for the real-space matter correlation function 
within the Zel'dovich approximation. It has the apparent form of a 
Gaussian convolution kernel, except for the fact that the matrix Ay 
is a function of q. Indeed, we see that £ <ZA) arises entirely from the 
scale-dependence of this Lagrangian correlator. 

The smoothing of the acoustic feature is often modeled as a 
convolution of the linear correlation function by a Gaussian kernel, 
with the smoothing scale estimated at lowest order by 2cr 2 . We can 
massage our expression into a similar form by noting that Ay can 
be written as the sum 

Aij(q) = By + Qjiq), (35) 

where By = 2cr 2 8ij is scale-independent. Then, by the same rea- 
soning as is used to show that the convolution of two Gaussians is 
a Gaussian, we can write 

where we have defined 

1 = J {2 jL,n e-^-^-»K (37) 



)3/2| C |l/2 



Eq. d36b is a proper Gaussian convolution, since the matrix B is in- 
dependent of q. The quantity x(q) ma y therefore be viewed as an 
analog of the linear correlation function. Indeed, the two are quite 
similar. These observations provide analytic jus t ificati on to con- 
ventional wisdom, first pointed out in iBharadwail ( fl99^) , that non- 
linear structure growth causes a Gaussian smearing of the cluster- 
ing signal. In our approach, this result is obtained at leading order, 
within the Zel'dovich approximation. 



3.2 Perturbative expansion for biased tracers 

Returning to the case of biased tracers, consider again Eq. l |27t . In 
the unbiased case the k integration in Eq. \20\ took the form of a 
Gaussian integral, which we carried out analytically. In the biased 
case, we can achieve the same thing if we first partially expand Eq. 
j27l l as 



K 



[l-A 1 A 2 ^ R -(A l +A 2 )Uik 



+ -A]A\f R + ^(A l +A 2 ) 2 U l Ujk,k j 
+ AMAi + Ai^RUih + 0(P 3 L )] 



(38) 



We may justify this choice of expansion by noting that both £n(q) 
and Ui(q) vanish in the large-scale limit |<jr| — » 00, while a 2 R and 
Ajj(q) approach non-zero values. In the language of perturbation 
theory, keeping these terms exponentiated therefore amounts to an 
non-perturbative resummation of the dominant large-scale contri- 
butions. 

iTassev & Zaldarriagal J2012al lbl) have recently emphasized the 
importance of not splitting the effects of bulk flows across orders 
in perturbation theory. The resummation described above has this 
proper ty, which is not shared by th e resummations used in LRT or 
RPT ( Crocce & Scoccimarroll2006h . 

To get from Eq. d38b to an expression for £x(r), we must in- 
tegrate K over Ai,A 2 , k, a nd q. The Ai and A 2 integrations may be 
evaluated via the identity (Matsubara 2008j3) 



r dA 

J lH 



F(A) (iA)" e~v* 



dS -^^/ FW \ 



V2tt. 



.ncr R 



dS" 



(39) 

where (^F M \ is the expectation value of the nth derivative of the 
Lagrangian bias function F(S) (see Appendix[A]for details). Appli- 
cation of this identity leads to 



L(q,k) 



/dA\ dA 2 - 
— T-j^ F(A l )F(A 2 )K(q,k,A u A 2 ) (40) 

1 



l+{F'?t;R + 2i{F')U l k l + -{F") 1 t;l 



- «F"> + {F') 1 )U,U j k,k J + 2i (F') (F")&Ufc 
+ 0{P])] . 



(41) 



The k integration reduces to a series of multi-variate Gaussian in- 
tegrals of the form 



d 3 k 
(2nf 



(42) 
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Appendix Icl reviews the relevant formulae. In the end we obtain 

M(r,q) ee [ i!L e *<*-'>L(g,fc) (43) 
J (2n) 3 



1 



(2w) 3 / 2 \A I " 2 



[l+<^"> 2 ^ 



1 



- 2<F'> U,g, + - (F") 2 f R - ((F") + (F'fWtUjGij 



where 



2{F'){F")% R U igi + 0(Pi)], 



g^A'iq-r), (i.-iA hj-gtgj. 



Our final expression for the correlation function is 



l+&(r)= d 3 qM(r,q). 



(44) 



(45) 



(46) 



The remaining integration over q must be performed numerically. 

Note we l l that, a lthough our calculation is very similar to that 
of Ma tsubarH fe008bl) . our result for £x(r) is not simply the Fourier 
transform of his Eq. (34) for P„bj(k) The difference lies in our 
choice of expansion in Eq. d38t . As discussed previously, the matrix 
Ajj(q) is the sum of a constant t erm 2o -2 ,<5,, and a scale-dependent 
remainder Cjj(q). In Matsubara (2008b) only the constant piece is 
exponentiated while the rest is expanded, i.e. 

^Mat = e-^ +A >le-^ kl 1 - ^k T Ck - AiA^ R - (A, + A 2 )U,k, + ■ ■ 

(47) 

Our approach may be seen as a partial resummation of the result of 
iMatsubad |2008b), and as such we expect it to be more accurate 
on small scales. 

Before we leave this section it is worth noting the manner in 
which the bias terms enter in Eq. d4 1 b . In particular note the term 
which goes as {F'){F") at the end of the 3 rd line and the (F") 2 
term at the end of the 2 nd line. For highly biased halos, assuming 
the peak-background split to compute the bias, {F") oc (F') 1 oc b 2 , 
so these terms can come in with (apparently) large powers of b, 
beyond the b 2 ter ms which one would naturally expect in a 2-point 
function (see also Reid & WhitefcoTll) . In our calculation using the 
Zel'dovich approximation and local Lagrangian bias we see these 
important contributions arise from 2 nd order bias. 



3.3 Redshift space 

Thus far we have concentrated on real space results, however the 
transition to redshift space is easily achieved. Recall that the dis- 
placement field in redshift space is given by 'V s = *P + H~ l (z ■ *P)|. 
In the Zel'dovich approximation *P oc D(t), so 



(ZD 



(48) 



Our previous derivation remains valid, we need only make the sub- 
stitutions 



U,^ U* = (S, j +fz,Zj)U ] , 



Aij 



(Stk + fZiZk)(Sji + fZjZi)A u . 



(49) 
(50) 



This slightly complicates the evaluation of the q integration in Eq. 
d46l l. in that we can no longer use azimuthal symmetry to reduce it 
to a 2-D integral. Nevertheless, the full 3-D integral is still feasible 
numerically, and the redshift space correlation function t; s x (s) may 
be easily calculated. 



3.4 Linear theory limit 

Standard Eulerian perturbation theory describes an expansion for 
the power spectrum of the form, 

P(k) = P {l \k) + P i2) (k) + ■■■ (51) 

where Z*'' is 0{P" L ). Unfortunately this expansion does not translate 
into a well-defined perturbative expansion for £(r), as the Fourier 
transform of P ( " } diverges for n > 1. Nevertheless, the linear the- 
ory correlation function is well-defined, and our approach should 
reproduce this limit when P L is small. We now show that this is 
indeed the case. 

In the Zerdovich approximation, the correlators Ay and {/, are 
given by linear integrals of P L , 

■01 - - -kkj 



C d 3 k r .. 



k 4 



-Pdk), 



(52) 



(k). 



The quantity M(r, q) defined in Eq. < !44t > is therefore ill-defined in 
the limit P L — > 0. To make our discussion precise, we replace the 
matrix A,j in this expression by 

B, J (q)=f3 2 6, J + eA, J (q), (53) 

where {S is a regularizing parameter that will eventually be set to 
zero, and e is a book-keeping parameter to help keep track of pow- 
ers of P L . Thus we write 



1+fxW 



lim I — 



\B\ 



e -4(/- ? ) T B-'(/-9) M + 



26<F'> l/ T B-y - q) + e{F') 2 £ L + 0(e 2 )] (54) 

(55) 



= \ + e^\r) + e 2 £\r) + --- . 

The linear contribution is then given by = d^ x /de\ e=0 . 
Using the identities 



d 

— detB = (detB)Tr 
ae 



ae 



as- 1 



we have 



d§x f d 3 q 

de J (2nfl 2 \B\ 



-B — B 

ae 



4(r~4) T B-'C--4) 



(56) 
(57) 



1/2 



i(r - 9 ) T B-'AB-'(r -q)-\ Tr(B _1 A) 
+ 2 0F'> U T E- > (r-q) + (F') 2 ^ L + 0(e)] 

~ J (27r)3/2/?3 E [2 P (r q) Mr 

- l -IT 2 Tr A + 2 (F')IT 2 U T (r - q) + (F'f ft 



Integrating by parts, and noting that 
1 



£S> (2/r) 3 / 2 /3 3 



(r-q)-llp l 



■ 6 D (r - q), 



we obtain 



1 d 2 A u , dU, , , 

o T-Sr-W ~UF') -^(r) + {F'ftdr). 

2 orjdrj orj 



We see immediately from Eq. \52\ that 



d 2 A i: 



-(q) = 2fe(ff), -^(q) = -Uq\ 
aqiOqj aq { 



(58) 



(59) 



(60) 



(61) 



(62) 
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Therefore the linear theory limit of our result is 

Q\r) = [\ + {F')ft L (f), (63) 
in agreement with standard perturbation theory. 



4 HIGHER ORDER 

We now repeat the derivation of the previous section, this time with 
the aim of extending our result to one order beyond the Zel'dovich 
approximation. Many of the technical details are relegated to ap- 
pendices. 

We pick up the track following Eq. J2U . prior to which we 
make no use of the Zel'dovich approximation. With the help of the 
multinomial theorem, the cumulant expansion of Eq. \2\\ in the 
general case can be written 

log K = E 7JT ( (Al6> + Xl&1 + k ' A)N )c (64) 

N=\ 

■ni+n+r 

= X -mWh ■ ■ ■ K (WA,, . . . A„.) . (65) 
^— ' m\n\r\ x /c 

m,n,r 

The cumulants (SfS^A^ . . . A;,\ are the key ingredients in our the- 
ory. In the following we refer to them generally as "Lagrangian cor- 
relators." As emphasized previously, they are functions of q only, 
so their tensor structure places severe restrictions on their func- 
tional form (see Appendix [Bj. Moreover, due to the properties of 
Gaussian random fields, a cumulant of order m + n + r must be at 
least of order m + n + r — 1 in the linear power spectrum P L (e.g. 
iBernardeau et al.l[2002h . An expansion in cumulant order therefore 
corresponds to a perturbative expansion in powers of Pt- 

For convenience, we assign different symbols to these La- 
grangian correlators based on their tensor rank r. For r = 0, since 
Sr is Gaussian, the only non-vanishing cumulants are 



(8\) c = (Sl) c = a 
For r = 1,2, and 3 we denote 



{SiS 2 ) c = &(«). 



W™ ee <<^A,AA>c 



(66) 



(67) 
(68) 
(69) 



Explicit expressions for these quantities may be found in Appendix 
|B1 Since they arise frequently, and to remain consistent with the 
previous section, we also adopt the shorthand 

Uj° -> U„ A™ -> Ay, and Wf° k -> W ijk . (70) 



In this notation, we evaluate Eq. J65t up to cumulants of order 



three, 
log AT 



(At + A 2 )U,k, - -W l]k k,kjk k - -(At + A 2 )Alfk i kj 
o z 

l -(A\ + A^Ufki - iA { A 2 Ul l k, + O(Pl). (71) 



We recover K by exponentiating. Of the eight terms in the above 
expression, only the first two have non-zero limits as |<jr| — » oo, and 
include O(Pt) contributions. As in the Zel'dovich case, we leave 



these two terms exponentiated while expanding the rest, thus 
K = e -¥A +A >l-i A iM^\ - A { A 2 % R - (A, + A 2 )Ujkj + ^A\A\f R 
+ l -(A { + AofU.Ufrkj + AiA 2 (Ai + A 2 )t R U,k, 
- l - WijMjh - l - + A^Ap.kj - l -(X\ + A 2 2 )Uf°k, 



U^Uj'k + 0(P|)]. 



(72) 



As in Section L3T21 we must now integrate with respect to A\ 
A 2 , k, and q. The analog of Eq. J4U for the one-loop case is 



1 + <F'> 2 fi, + 2i (F') Uih + - (F"f f R 



- «F"> + (F'rW.Ujk.kj + 2i (F') {F'^^Uih 

- Wijkhkjh - (F) Ajfkikj + i {F") Uj% 



+ i(F') 2 Uj l k, + 0(Pl) 



(73) 



Analogous to Eq. J44t . the k integration gives (see Appendix let 
1 



M 



(27r) 3 / 2 |A| 1 / 2 



-4(i-?) J A-'(r- 4 ) 



1+W & 



2 (F') U igi + I (F") 2 & - [<F"> + (F') 2 ](JjU jGjj 



- 2{F') {F")^ R U,g, + -U.J.. - (F')A}fGij 

- (F'^Ufg.-iF'fuj'gi + OiPl) . 

where g, and Gy are defined in Eq. j45\ , and 

r, 7 * = (A -1 ) y gt + (A- l ) kigj + (A- l )j k gi - g igj g k , 



(74) 



(75) 



Our final expression for the real-space correlation function £x( r ) is 
given once again by Eq. l !46t . with M(r, q) given by Eq. d74b up to 
0(Pl). The redshift -space correlation function §(s) is obtained by 
replacing the real-space Lagrangian correlators by their redshift- 
space counterparts. 



5 RESULTS 

Having presented the formalism and rationale behind our resumma- 
tion, we now compare the results of our "convolution Lagrangian 
Perturba tion Theory" (CLP T) to linear theory and to the earlier 
work of lMatsubaral d2008allbl) . This is the most natural compari- 
son, since our work is largely an extension of LRT and a partial 
resummation of that formalism. 

Fig. [U shows the (monopole) matter correlation function 
in real- and redshift-space. The solid line shows linear theory 
while the dashed and dotted lines show our CLPT and Matsub- 
ara's LRT respect i vely. In redshift-space we have used the for- 
malism of Kaiser ( 19§3) as our "linear t heory". The po i nts ar e 
from the N-body simulations described in iReid & White] d201lh : 
IWhiteetal.l ( [201 ]]). Throughout this paper we compare exclusively 
with z = 0.55 simulation outputs. Note that linear theory provides 
a poor approximation near the peak of the correlation function (at 
100/r'Mpc) in both real- and redshift-space, as is well known and 
we have discussed previously. On large scales LRT and CLPT are 
nearly indistinguishable, as expected. However on smaller scales 
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Figure 1. (Top) The real-space, matter correlation function, ^(r), from lin- 
ear theory (solid), LRT (dotted) and CLPT (dashed) compared to N-body 
simulations (squares) at z = 0.55. In order to plot the results with a lin- 
ear y-axis we have multiplied f by r 2 , which removes much of the trend 
from r ^ - lOOMpc. LRT and CLPT agree very well on large scales 
(the lines can barely be distinguished) and agree well with the N-body 
results. LRT overshoots the N-body results below r 20/i~'Mpc while 
CLPT tracks the N-body results to much smaller scales. Linear theory over- 
shoots at r =s 20/r'Mpc and at r ^ lOO/r'Mpc. (Bottom) The redshift- 
space, monopole, matter correlation function, foM, from linear theory 
(solid), LRT (dotted) and CLPT (dashed) compared to N-body simulations 
(squares). The qualitative behavior is as for f (r). 



the resummation inherent in our approach allows CLPT to track 
the N-body results to smaller scales than LRT. 

The comparison with the quadrupole and hexadecapole mo- 
ments of the redshift-space correlation function is very similar 
(Fig. [2}. Both LRT and CLPT provide a better fit than linear theory 
to the quadrupole and hexadecapole moments at large scales, but all 
theories depart from the N-body results at larger scales than for the 
monopole. The level of agreement is worse for the hexadecapole, 
but that moment is also quite small. 

Fig. [3] compares the theories for biased tracers, in this case 
for halos in the range 12.8 < lgM;, < 13.1 at z - 0.55 though 
other results are qualitatively similar (see Fig. |4jl. The situation is 
similar to that for the matter: linear theory provides a poor approx- 
imation at large scales, missing the smearing of the acoustic peak 
due to the motion of material. LRT tends to overshoot the N-body 
results at small scales, while CLPT provides a good match down 
to (9(10/i~'Mpc). Note that we considered two distinct sets of bi- 



Figure 2. The redshift-space, quadrupole and hexadecapole, matter corre- 
lation functions, and £i(.s), from linear theory (solid), LRT (dotted) 
and CLPT (dashed) compared to N-body simulations (squares) at z = 0.55. 
For the quadrupole LRT and CLPT agree very well on large scales (and 
agree well with the N-body results) but LRT departs from the N-body re- 
sults at much larger scales. For the hexadecapole the disagreement between 
N-body, CLPT, LRT and linear theory breaks down at larger scales than for 
the quadrupole. 

asing parameters. In Figs.[5]and[5]we allowed the "renormalized" 
bias parameters (F'> and (F") to be adjusted independently, while 
in Fig.[H we relate d the two using the peak-background split, as in 
iMatsubanJ d2008bll3> . 

Finally we compare the monopole and quadrupole moments 
of the redshift space correlation function of halos to the predictions 
of CLPT in Fig. [5] The prediction of the monopole moment is in 
relatively good agreement with the N-body measurements, though 
the level of agreement at ~ 20/r'Mpc is clearly not as good as 
it was with the matter. The prediction for the quadrupole is much 
worse than it was for the matter. 

On large scales the prediction for the quadrupole is dominated 
by the same terms as the matter and the term scaling as {F'). The 
CLPT prediction does not have as much power on small scales 
as the N-body results, which have more small-scale power com- 
pared to the large-scale power than was the case for the matter. 
The shortfall in power is shared by the terms which survive when 
(F') = and by the terms which scale as (F'). The failure of 
our model to match the quadrupole moment on small and inter- 
mediate scales may be due to our assumption of local Lagrangian 
bias. While this approximation has received some support from 
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Figure 3. The real-space, correlation function for halos with 12.8 < IgA//, < 
13.1 computed in linear theory (solid), LRT (dotted) and CLPT (dashed) 
compared to N-body simulations (squares) at z = 0.55. In this plot we al- 
lowed (F 1 ) and (F") to vary independently to obtain the best agreement 
with the N-body results. 
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Figure 4. The real-space, correlation function for halos in three mass bins 
computed in linear theory (solid), LRT (dotted) and CLPT (dashed) com- 
pared to N-body simulations (squares) at z = 0.55 for three different mass 
ranges each a factor of two in width: from bottom to top 12.2 < lgM/, < 
12.5, 12.8 < lgM A < 13.1 and 13.1 < lgM,, < 13.4 with masses in h" l M Q . 
In this plot we enforced the peak-background split relation to determine 
(F") in terms of the best fit (F'), i.e. the theory has only one free parame- 
ter. 



N-body simulations dRoth & PorcianH 1201 lUBaldauf et all |2012| ; 



IChan, Scoccimarro. & Sheth 20121; IWang & Szalavll2012^ we also 
expect that terms involving e.g. the tidal tensor, can become impor- 
tant for high mass halos JSheth, Chan & Scoccimarroll2012h . Such 
terms are naturally quadrupolar in nature and may affect the pre- 
dictions. 



6 DISCUSSION AND CONCLUSIONS 

We have presented a new formulation of Lagrangian perturba- 
tion theory which allows accurate predictions of the low-multipole, 
real- and redshift-space correlation functions of the mass field and 
dark matter halos. Our formulation, which we refer to as "convo- 
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Figure 5. The redshift-space, monopole and quadrupole, correlation func- 
tions for halos computed in linear theory (solid), LRT (dotted) and CLPT 
(dashed) compared to N-body simulations (squares) at z = 0.55. 



lution Lagrangian perturbation theory" or CLPT involves a non- 
perturbative resummation and in deed can be view ed as a partial 
resummation of the formalism of Matsubara (2008a b) upon which 
we have relied heavily. 

A key difference between CLPT and LRT or RPT is that 
we naturally recover the Zel'dovich approximation as the low- 
est order of our expansion for the matter correlation function. 
iTassev & Zaldarriagal d2012al) have recently emphasized the impor- 
tance of not splitting the effects of bulk Hows across orders in per- 
turbation theory, and we find that CLPT (which does not make such 
a split) does indeed provide better agreement with N-body results 
at small scales than LRT (which does). 

CLPT works best for the real-space clustering of the matter 
and halos and for the monopole of the redshift-space correlation 
functions. While the N-body results for the quadrupole and hex- 
adecapole moments of the redshift-space correlation function for 
the matter is relatively well reproduced by CLPT, those moments 
for the halo correlation function differ significantly from the CLPT 
prediction. We suspect that this difference is due to a limitation in 
our bias prescription, in particular that our assumption of local La- 
grangian bias for halos is not sufficiently accurate. Further work 
along these lines is clearly warranted. 

One possible extension of this work is to use the real-space 
correlatio n function from CLP T in the Gaussian streaming model 
ansatz of iReid & White] d201ll) with vn and 0" 12 terms calibrated 
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from N-body simulations or computed within the context of LPT. 
These terms can be computed in our formalism by generalizing our 
function K (Eq. I2U to include a A contribution and taking func- 
tional derivatives of K. We leave this for future work. 

Finally, we note that our work may be relevant for efforts 
to model the bispectr um within the Lagrangian framework (e.g., 
iRampf & Won32012h . 
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APPENDIX A: BIAS 

As in lMatsubarj |2008b). we note the identity 

— F(X)e-^ Al ' rl K{iX)" = (F M ) , (Al) 

where \F^ n A is the expectation value of the nth derivative of the 
Lagrangian bias function, what are referred to as "renormalized" 
bias coefficients c„ in Matsubara 2011. The mapping K — > L is 
therefore achieved by replacing 

(Ai+A 2 )-»-2/<F>, (A2) 

A^i^-iF') 2 , (A3) 

4^-><F"> 2 , (A4) 

(A l+ A 2 ) 2 ^-2[(F") + (F') 2 ], (A5) 

A l A 2 (A l +A 2 )^2i{F')(F"), (A6) 

x \ + A l -» -2(F"). (A7) 

APPENDIX B: LAGRANGIAN CORRELATORS 

In this appendix we collect the relevant facts and formulas concern- 
ing Lagrangian correlators that we need for our one-loop theory. 
The correlators are defined by 

C™.„.(?) = (ST^At, -Afr) e , (Bl) 
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where 6 { = 8 L (q{), 8 2 = 6 L (q 2 ), and A, = T ; (^ 2 ) - The 
subscripted c refers to a connected moment; since these fields have 
zero mean, the connected moments coincide with normal expecta- 
tion values for orders m + n + r < 3. 



Bl Index structure 

By translational symmetry, a Lagrangian correlator can only de- 
pend on the Lagrangian separation q = q 2 - q l . This imposes 
strong constraints on its index structure. We classify a correlator 
by its tensor rank, i.e. by the number of vector indices it carries. In 
the following we let t/,, Ay, and Wjj k denote generic correlators of 
ranks 1, 2, and 3, respectively. 

Rank- 1 correlators must be of the form 



U t (q) = U(q) q, 



(B2) 



for some scalar function U(q), since (trivially) the only vector 
quantity that can be formed from the vector q is proportional to 
q. Rank-2 correlators must involve only rotationally invariant rank- 
2 tensors that can be formed from the vector q, i.e. Sjj or q^j. Thus 
their general form is 

A,j(q) = X(q) Sij + Y(q) qi qj. (B3) 

Likewise, rank-3 correlators are of the form 

Wijk(q) = Vi(q) qiS jk + V 2 (q) qj8 ki + V 3 (q) q k 8ij + T(q) Mj9k- (B4) 

We remind the reader that we adopt the shorthand 

£/, 10 ->£/„ A™->A, P and -» W i;h (B5) 

since these combinations arise frequently. 

In general, correlators of even rank are even functions of q, 
while those of odd rank are odd. This implies that the correlator 
C™" . is symmetric in the indices m and n, as the following chain 
of equalities shows: 

C™,,(?) = (^A„...A, r ) c 

= {6L(q l )"'6dq 2 )"W,M 2 )- x V,M^ • • • 

= {-l) r {SL^TSdqxT^tMO - ^(ffa)] • • • (B6) 
= (-l) r C™".„.(-9) 

We can solve for the coefficients in these expansions by con- 
tracting against tensors and solving the resulting simultaneous 
equations, e.g. for the components of W,^ : 



3V l + V 2 + V 3 + T= W, jkqi 6j k , 
V 1 + 3V 2 + V 3 + T= W, jk QjSu, 
V [ + V 2 + 3V 3 + T=W, jk q k S, J , 
V l + V 2 + V 3 + T = W, jk q,q iqk . 

B2 Perturbative orders 

The LPT expansion of the field A has the form 



(B7) 



where A <a) involves a factors of the linear density field <5 L . The cor- 
relators C™. , may therefore be expanded as 



(B9) 



a\ = l a,=\ 



where Cy 



' r) = (dfSlA!? ■■■A < ? r) \ . Since 6 L is Gaussian, 
many of these terms vanish. Here we display the breakdown for 
each of the quantities introduced in Section[4] up to order O(Pj): 



Ui 


= f/ a) + t7 (3) + ... ; 


(BIO) 




= A oi) +A a2) +A (i3) +A ai) + ... j 


(Bll) 


w, Jk 


ijk ijk ijk ' 


(B12) 


u? 


= t/ 20<2) + ..., 


(B13) 


up 


= y, I1<2) + ---, 


(B14) 


A- 




(B15) 



B3 Scalar components 

Given the index structure described in the previous subsection, 
evaluating the Lagrangian correlators reduces to computing a set 
of s calar functions of q . In order to maintain notational consistency 
with lMatsubari (2008b) we make use of his definitions of Q and R. 



° L {kr)R n (r) 



R„(k) = ■^■Pdk) J drP 
and 

Qn(k) = ^J drPdkr) j \lxP L (k^)Qn(r,x) , (B 
where y(r, x) = 1 + r 2 - 2rx and the Q„ are given by 



(B16) 



17) 



Qa 

Q 3 



2 (l-* 2 ) 2 



Qi 



(1 - ^rxil - rx) 



and 



T T 

x 2 {\ - rx) 2 1 - x 2 

~f ' ^ 4 ~ ~f ' 

rx(l-x 2 ) (1- 3rx)(l - x 2 ) 

<Js - ■ > <J6 ~ > 

^(1 - rx) r 2 (l - J) 

Ui - . <J$ = . 

y y 

rx(l - rx) , 

Qg = , Qw = 1 - JT, 

7 y ■> 

Qn=x~, Qn = rx, Q [3 =r 2 



i r 2 (l-x 2 ) 2 



R l( r) = £dxg 



r2 - Irx 



- p+ 1 , (1 x 2 )rx(l - rx) 

Mr) = J_, dx — — - — 

J 1 1 + r — 2rx 

In the follo wing, equat i on refe rences prefaced with "M" indicate 
equations in lMatsubari] d2008bh . 

The expression for Ay = A?? is derived in detail below. The 
other components we need are 



A ( " + A (2 ' + A 



(3) 



(B8) 



= <*iA,Aj> e 

= Xi (q)Sij + Y 10 (q)qiqj 



(B18) 
(B19) 
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with 



V\ U2 \ q ) = 

v< 112, (^ 

S (U2) (^ 



t/ <3) (9) 



<(<?) 



^ J dkeP L (k)j {kq) 

iff ? [M ' + "' + ft + 2|a 

1 f°° <ft / 3\ 

2^ Jo ~ rvj x 

[2Ri+4tf 2 + g 1 +2fi 2 ]./3(fy) 

-i- f dkk{-l)P L {k)h(kq) 
2n- Jo 

2^1 
^X 



|-^lfiByi(*?) 



>; f) ^--|[J?i +R2]jl(kq) 



-3[3R, + 4R 2+ 2Q 5 ] J ^\ 



2n 2 J 



dk \-—\[3Ri + 4R 2 + 2Q 5 ]x 



jo(kq) - 3 



Mkq) 



(B20) 
(B21) 
(B22) 
(B23) 

(B24) 
(B25) 

(B26) 

(B27) 

(B28) 

(B29) 

(B30) 



where the arguments of the R„ and Q„ terms are k and have been 
omitted for brevity. The remaining equations, for X^ n \ X^ 22 \ Z (13) , 
F (22) , y (13) are presented and derived in the next section, 



B4 Example 

We provide here an example of how to obtain the formulae of the 
previous subsection. We focus on Ay = (A;Ay) c , since this is the 
most important of the Lagrangian correlators in our theory. 
By the definition of A, 

A, =>F,(? 2 ) -¥,(?,) = J |^(e' >42 -e iMl ) (B31) 
and therefore 

x (v.ip^M)}^ (B32) 
From Eq. (M.A9), the Fourier space 2-point function here is 



J (2^)3 (2n?\ )\ e > 



{ x Vi(P l ) x i'j(P2)) c = -(2n?Sl( Pl + pdCfa). (B33) 

The quantity Cij(k) here has contributions at both tree and l-loop 
level, 



cf\k) = c^ k ) = -^-Rm. 



(M.A52) 
(M.A53) 
(M.A54) 



These terms are all of the form Cy = -{kjkjlk^)a{k) for scalar a(k), 
as is guaranteed by rotational symmetry. With the substitution of 
Eq. dB33t into Eq. ( |B32t . we have 

Aii = I ^ 2 - e ' k ' q - e " 9 ) kj w a{kr (B34) 

Contracting this quantity first by <5y and then by q,qj, we obtain the 
system of equations 



A iA . = 3X + r = JJ|( 
A, J q i q J =X + Y = |^(2-^-^) 



2 - e ,k " - e-' k ") —a(k), (B35) 
(k ■ qf 



k 2 



-a(k). (B36) 



Letting fi = k • q and using the Bessel function identities in Ap- 
pendix[C]\ve may perform the angular integrations, 

+ Y = JL j k 2 dk l - jf dfi (2 - e< k <>» - e*») ±a(k) 
= iX ^t 2-2 ^*?)]^*). ( B3? ) 



2^X 



(ft 



- - 2] (kq) + 4— 

3 kq 



a(k), (B38) 



from which we obtain 

2tt- Jo 
2tt 2 Jo 



2 2 Ji(%) 

3 fc 9 



-2j Q (kq) + 6 



fa/ 



(B39) 
(B40) 



Explicitly, up to l-loop order, the contributions to X(q) and Y(q) 
arc 



dkP L (k) 



(k) 



2n 2 X 

X{ " Kq)= 2h[ dk 2l R ^ 
dkP L (k) 



2 j\ikq) 

3 kq 



2 _ j\(kq) 

3 kq 

2 Mkg)' 

3 kq 



Y m (q) : 



2^ 2 Jo 



-2jo(kq) + 6 



Mkq) 



Y<22>iq)= ^r dk h Qi(k) 



Y™(q) 



1 f°° 5 

5?JL *2i* 



(A) 



-2jo(kq) + 6 
-2j (kq) + 6 



kq 

Mkg) 



kq 
Mkq) 



kq 



(B41) 
(B42) 
(B43) 
(B44) 
(B45) 
(B46) 



Note that each of these quantities approaches as q — > 0. 



APPENDIX C: REFERENCE FORMULAE 
CI Gaussian integrals 

In our theory we make use of the basic Gaussian integral 



Q(b) 



r d 3 k 

J (2^F 



k -^k T Ak+ibk _ j „-U r A-'ft 



(2n-) 3 / 2 |A| 1 / 2 



(CI) 
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where \A\ denotes the determinant of the 3x3 matrix A. By succes- 
sive applications of the operator -id/dbj, we also have 

f -^G(k)k, = i{A- l b) t Q(b) , (C2) 

= [(A-'^-fA-'iMA-'^jeW, (C3) 

= /[(A-'MA- 1 *)* + (A- 1 )^- 1 *), 
+ (A- 1 ) yi (A- 1 A) i 

- (A-^MA-^J/A" 1 *)^*). (C4) 



r d 3 k 

J (2/r) 3 J 
) 3 



P rf 3 / 



where we have written 

G(k) = e -ifA* + ».* 
for notational compactness. 



(C5) 



C2 Spherical Bessel functions 

In performing the integrals in the previous sections we have found 
the following spherical Bessel function identities and integrals to 
be useful: 



±J><* 

2"£^ e "< 



(2n + 1) 

X 


(C6) 


{In + l)—j n (x) 
ax 


(C7) 


jo(x) 


(C8) 


iji(x) 


(C9) 


1 2 


(CIO) 


;oW-2 

X 


(Cll) 


3 2 
« jjiw - 5 Aw 


(CI 2) 



